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ABSTRACT 

We present high resolution mid-infrared (mid-IR; 11.7 and 17.65 fim) maps of 
the central 400 pc region of the starburst galaxy M82. Seven star forming clusters 
are identified which together provide ~ 15% of the total mid-IR luminosity of 
the galaxy. Combining the mid-IR data with thermal radio measurements and 
near- and mid-IR line emission, we find that these young stellar clusters have 
inferred masses and sizes comparable to globular clusters. At least 20% of the 
star formation in M82 is found to occur in super-star clusters. 

Subject headings: galaxies: individual(M82) — galaxies: starburst — galaxies: 
star clusters — infrared: galaxies 



1. Introduction 

M82 (NGC 3034) is often considered the archetypical starburst galaxy since it has a 
derived star formation rate (~10 Mq yr~^; O'Connell & Mangano 1978) that would deplete 
the observed molecular gas in < 10® yrs (i.e. short on a Hubble timescale; Lord et al. 1996). 
The starburst phenomena traces recent star formation since it is the massive stars (>8 Mq), 
which have short lives, that dominate the energetic output of the host galaxy. The nuclear 
starburst of M82 dominates the infrared (IR) luminosity of the galaxy - essentially all of the 
galaxy's ~ 3.6 x 10^° Lq^ comes from the central kpc (Telesco & Harper 1980). 

Optical and near-IR imaging of M82 has revealed numerous super-star clusters in its 
active star-forming nucleus. Ground-based optical imaging of the central region detected 
eight young knots distributed throughout the region (O'Connell &; Mangano 1978). Further 



^Telesco & Harper (1980) published Ljr = L(l-300/im) = 3.0 x 10^° Lq assuming a distance to M82 of 
3.3 Mpc. We have updated Lju for the distance of 3.6 Mpc (Freedman et al. 1994; Sakai & Madore 1999) 
used in this paper. 
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study of these knots broke them into smaller star clusters (half hght diameters ~3-4 pc) 
and suggested cluster dynamical masses 10^-10^ Mq (Smith & Gallagher 2001; de Grijs, 
O'Connell, & Gallagher 2001; McCrady, Gilbert, & Graham 2003). 

M82's current starburst is thought to have been triggered by its interaction with M81 
~10^ yrs ago (Cottrell 1977; Achtermann & Lacy 1995). It has been argued that the con- 
ditions resulting from interactions and mergers of galaxies are favorable for globular cluster 
formation (eg. Taniguchi. Trentham, & Ikeuchi 1999 and ref. therein). UV studies of global 
properties of starbursting galaxies have shown that as much as 20% of the light is produced 
in the luminous knots (Mcurcr ct al. 1995; Zcpf ct al. 1999), suggesting a high efficiency of 
cluster formation in starbursts. Understanding the star formation occurring in M82's nu- 
clear region can provide insight into both the general process of star formation in starburst 
environments and also the process and efficiency of forming super-star clusters. 

The nearly edge-on geometry of M82 combined with heavy optical extinction [Ay ~ 
5-25 mag; Lester et al. 1990; Telesco et al. 1991; Larkin et al. 1994; Satyapal et al. 1995) has 
made the galaxy's central 400 pc difficult to study. Even in the near-IR where A2.2//m ~ 
1/lOAy (Rieke & Lebofsky 1985), it is difficult to directly measure the deepest star forming 
regions. The mid-IR region of the spectrum is essential for probing deep into the central 
regions of M82 and revealing the details of the intense star formation occurring in the central 
regions. Previously, the highest resolution mid-IR map of the central region of M82 was a 
12.4 yum map by Telesco & Gczari (1992), which had a resolution ~1.1". Mid-IR maps of 
M82 with lower resolution have also been published by Rickc ct al. (1980), Telesco, Dcchcr, 
& Joy (1989), Dietz et al. (1989), Telesco et al. (1991) and Fdrster-Schreiber et al. (2003). In 
this paper, we improve upon previous observations by presenting higher resolution mid-IR 
maps of M82. We discuss evidence suggesting the sources in the maps are young counterparts 
to globular clusters and estimate the efficiency of super-star cluster formation in M82. 



2. Observations 

On 2003 April 23, we imaged M82 at 3.5 (L-band), 11.7, and 17.65 /xm with the Long 
Wavelength Spectrograph (LWS; Jones & Puetter 1993), a facihty instrument at the W.M. 
Keck Observatory. LWS uses a 128 x 128 pixel Boeing Si:As detector and has a plate scale 
of 0'.'08 pixel"^, resulting in a 10'.'2 x 10'/2 field of view. We used the "chop- nod" mode of 
observing, with a chop throw of 15" north. The bad pixels in the images have been smoothed 
over, and a mask has been applied to remove the portion of the chip not illuminated by the 
source. The seeing varied during the course of the observations so the resolution in individual 
frames ranged from 0.4-1.0" at 11.7 /xm and from 0.5-0.7" at 17.65 /xm. At each of seven 
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pointings across M82's nuclear region, we imaged at all three wavelengths before moving 
to the next pointing. Images from the seven pointings were mosaiced by centroiding on 
the bright sources in each frame. Assigning coordinates to the field was accomplished by 
aligning 2MASS sources with the centroids of the two bright sources in the L-band frames. 
We estimate that our positions are accurate to ~0'.'5. The star /i UMa was used as the 
primary standard for flux density and point-spread function calibration, and a Her and rj 
Sgr were used to estimate a cahbration error of ~ 20%. 

3. Chciracterization of mid-IR sources 

Our mid-IR maps, presented as Figures 1 and 2, contain several bright, resolved sources 
as well as diffuse emission connecting the brighter sources. The mid-IR sources, labeled 
A-G from west to east, are denoted by black circles (with radii O'.'S, corresponding to the 
positional error). For reference, the 2 /im peak is marked in the Figures as a yellow cross 
(Dietz et al. 1986) and the dynamical center measured by Lester et al. 1990 is marked with 
a yellow circle. The sharp edges visible in the maps (eg. southeast of source B in the 11.7 
fim map) are artifacts that resulted from the mosaicing of images with differing thermal 
backgrounds and do not affect the results of this paper. 

The earher, lower resolution 12.4 iim map by Telesco & Gezari (1992) contained two 
bright regions, one to either side of the galaxy's center. The overall structure in our higher 
resolution maps is comparable, though we identify individual resolved sources within the 
Telesco & Gezari sources. To compare their flux densities with those previously published, 
Telesco & Gezari reported a flux of 6.5 ± 0.7 Jy at 12.4 fim for the 4"x 4" region containing 
our sources C and D. We measure a flux for the same region of 5.3 ± 1.3 Jy at 11.7 /^m, in 
agreement with the previous observation. 

The flux densities and half intensity major and minor axes measured from the maps 
for the mid-IR sources are listed in Table 1. For sources covered by more than one image, 
photometry was performed on each image individually and the average of the measurements 
is presented here. The standard deviation of the measurements due to background variations 
was < 15% and, combined quadratically with the calibration error, results in an overall flux 
density error of < 25% in our measurements. 

Using our 11.7 /xm and 17.65 /^.m flux densities, we estimate the color temperature of the 
mid-IR emitting dust to be in the range 150-270 K, assuming the dust particles radiate as 
blackbodies. The mid-IR luminosities (Lmir = Li2-i8/im) of the sources, estimated by fltting 
the measured flux densities to a blackbody with the color temperatures for each source, are 
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all between 0.2-6 x 10^ and sum to 2.4 x 10^ L©. Lmir for the entire galaxy estimated 
using the same procedure with the uncorrected IRAS 12 and 25 flux densities of 53 Jy 
and 274 Jy, respectively, gives 1.7 x 10^° Lq, thus the seven mid-IR sources contribute ~15% 
of M82's total Lmir- 

Using radio images at flve frequencies, Allen (1999) created spectrally decomposed im- 
ages of thermal (free-free) and non-thermal (synchrotron) emission. We find that not only 
do all the centers of the thermal H II regions (marked in Figures 1 & 2 with red squares) 
match the centers of the mid-IR sources reasonably well, the diffuse structure in the mid-IR 
maps follows closely the structure in the thermal free-free map. Golla, Allen, & Kronberg 
(1996) also noted the correspondence of the diffuse 1.5 & 22 GHz emission with the Telesco 
& Gezari (1992) 12.4 fim emission. This supports the hypothesis that the mid-IR sources 
are heavily obscured H II regions. 

The [Ne II] (12.8 /xm) map published by Achtermann & Lacy (1995), with a resolution of 
~1", also correlates well with our mid-IR maps (see Figure 3). It should be noted that while 
the peak of the [Ne II] hue hes within our 11.7 //m filter bandpass (AA = 2.4 /xm), based on 
the 5-16.5 fjLm spectra presented in Forster-Schreiber et al. (2003), the [Ne II] line contributes 
<15% of the flux in our bandpass. We interpret the correlation between the mid-IR and the 
[Ne II] emission as confirmation that the source of the mid-IR emission is dust heated by 
UV from young stars which ionize Ne I. Further evidence to support this hypothesis comes 
from the Br7 map from Larkin et al. (1994). Though not covering the entire mid-IR field, 
the Br7 emission observed is also well correlated with the mid-IR emission. Both emission 
line maps agree with the mid-IR map in the apparent lack of emission toward the dynamical 
center of M82. 

Non-thermal radio sources from McDonald et al. (2002) and Allen (1999) are shown in 
Figures 1 & 2 as magenta crosses; these are assumed to be supernova remnants since they 
have inverted radio spectra. We find no correlation between the radio supernova remnants 
and the mid-IR sources, but notice that most of the supernova remnants follow the outer 
edge of the mid-IR emission at a flux levels <0.3 Jy arcsec"^ at 11.7 iim and <0.7 Jy arcsec"^ 
at 17.65 /im. 

Following Beck, Turner, & Gorjian (2001), we assume a reference OB star luminosity 
of 2.5 X 10^ L0 for an 07 star (Vacca, Garmany, & ShuU 1996) to calculate the OB star 
content of each mid-IR source from its Lmir- These values are listed in Table 1. For the 
brightest regions (B, C, and G) we find that ~2000-2500 07 stars are required to generate 
the mid-IR luminosity. Forster-Schreiber et al. (2001) derived comparable numbers of OB 
stars for regions near sources B and C (3200 and 4400 07 stars, respectively) using He I/Br7 
line ratios to estimate the Lyman continuum fiux from their regions. 
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Since we find that the [Ne II] emission correlates spatially with the mid-IR emission, 
wc use [Ne II] channel maps from Achtermann & Lacy (1995), which have a resolution of 
16 km s~^, to estimate velocity dispersions {a^^i) in the ionized gas of ~15-30 km s~^ at the 
positions of the mid-IR sources. Estimates of the velocity dispersions from the Br7 (Larkin et 
al. 1994) and ^^CO (Neininger et al. 1998) channel maps, each with lower velocity resolution, 
give similar results. A caveat to using ionized emission to estimate a velocity dispersion is 
that the ionized gas may well be accelerated by shocks, and therefore its velocity should be 
considered an upper limit to the stellar velocity. 

From the range of estimated velocity dispersions and assuming virializcd systems, we 
use M = rja^^irh/G, to calculate a range for the total mass in each these systems. In this 
formulation, G is the gravitational constant, Vh is the projected half intensity radius defined 
to be the geometric mean of the semi-major and semi-minor axes, and 77 = 10 (Smith & 
Gallagher 2001). The mass ranges for each mid-IR source are hsted in Table 1 and should 
be considered upper limits since H II regions are generally found to be freely expanding 
systems. For the largest sources, the mass range we find is 6-25 x 10^ M©. In every case, 
the range contains the system mass found by extrapolating from the number of 07 stars 
using a Salpeter mass function. Comparing the sizes and masses of the star clusters forming 
in M82 to three of the largest Galactic globular clusters, which have radii of 5-10 pc and 
masses 1-2 x 10^ M© (see Table 1), we conclude that the star clusters forming in the nuclear 
region of M82 are young analogs to globular clusters. 

The existence of present day globular cluster-sized knots of star formation is not unique 
to M82. Indeed, super-star clusters containing quantities of OB stars similar to those we 
find in M82 have been observed in a number of galaxies (e.g. NCG 5253, Turner et al. 2003; 
He 2-10, Johnson & Kobulnicky 2003; NGC 4038/9, Mengel et al. 2002). Additionally, two 
star forming regions in the Milky Way, the Arches cluster near the Galactic center and the 
Cygnus 0B2 association, are estimated to each weigh in at 6 x 10^ Mq (Serabyn, Shupe, & 
Figer 1998; Knodlseder 2000), comparable to a small globular cluster (Mandushev, Staneva, 
& Spasova 1991; Pryor & Mcylan 1993). However, the Arches cluster and Cygnus 0B2 each 
contain around 100 O stars, which contribute together only ~3% of the total Galactic O 
stars (Terzian 1974), and yet far outnumber all other Galactic star forming regions in their 
O star content. 

In M82, the seven mid-IR sources in the nuclear region together contribute ~ 15% 
of the total Lmir of the galaxy. The [Ne II] map suggests that there may be several H 
II regions outside our mid-IR field which may contribute up to an additional ~5% to the 
total LiR. Assuming all the mid-IR luminosity in M82 comes from star formation (Telesco 
1988), it follows that >20% of M82's star formation is in the form of super-star clusters - in 
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contrast to the mere 3% in the Milky Way. This may be an important feature of starbursts 
in general; not only do they provide an environment suitable for forming globular clusters, 
but the super-star cluster formation efficiency in starbursts is >20%. 

4. Conclusions 

This paper presents mid-IR (11.7 and 17.65 iim) maps with ~0.5" resolution of the 
central 400 pc of the star bursting galaxy M82. We find 7 resolved sources in this region of 
M82 with luminosities summing to 15% of the total IRAS flux of the entire galaxy. The 
mid-IR maps exhibit features comparable to those found in maps of [Nc II] emission, Br7 
emission, and thermal free-free emisson. We present evidence implying that the mid-IR 
sources arc giant H II regions in which globular cluster-sized star clusters arc forming. Our 
data imply that >20% of the star formation in M82 is occurring in super-star clusters. 
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Table 1. Properties of Proto-Globular Clusters in M82 and Galactic Globular Clusters 



Source 


F^(11.7 /xm) 


F^(17.65 /xm) 


Size'^ 


Tc 




N(07) 


Mass^ 




(Jy) 


(Jy) 


(pc) 


(K) 


(108 Le) 




(106 Mo) 


A 


0.16 


0.15 


10 X 10 


270 


0.20 


80 


0.3^ 


B 


3.17 


5.39 


30 X 21 


195 


5.9 


2400 


6-25 


C 


2.38 


5.95 


26 X 19 


165 


5.5 


2200 


6-22 


D 


1.23 


3.31 


21 X 12 


160 


3.3 


1300 


4-16 


E 


0.67 


1.74 


13 X 13 


160 


1.7 


700 


3-13 


F 


0.77 


2.24 


28 X 8 


155 


2.3 


900 


4-15 


G 


1.56 


4.96 


25 X 16 


150 


4.9 


2000 


5-20 


cu Cen 






9^ 








2^ 


NGC 6273 






3^ 








1^ 


NGC 104 






5^ 








1^ 



'^Sizes are full width at half intensity major and minor axes of elliptical sources, as observed, 
and include a ~0.5" (9 pc) PSF. 

'^Upper limit to mass of systems 

"^Source A is not identifiable in the [Ne II] channel maps. Stellar mass extrapolated from the 
number of 07 stars using a Salpeter initial mass function 

'^Sizes are averages of values in van den Bergh, Morbey, & Pazder (1991) and Mandushev, 
Staneva, & Spasova (1991) 

^From Mandushev, Staneva, & Spasova (1991) 
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Fig. 1. — Central region of M82 at 11.7 /im smoothed with a Of! 4 boxcar function. The seven 
mid-IR sources are labeled with white letters and their positional errors are represented by 
O'.'S radius black circles. The 2 /xm peak (Dietz ct al. 1986) is labeled with a yellow cross and 
M82's dynamical center as determined by Lester ct al. (1990) is marked by a yellow circle. 
Magenta crosses mark positions for non-thermal radio sources (McDonald et al. 2002; Allen 
1999) and red squares mark position of H II regions (Allen 1999). 

Fig. 2. — Central region of M82 at 17.65 //m smoothed with a 0'/4 boxcar function. Symbols 
are the same as in Figure 1. 

Fig. 3. — Smoothed 11.7 iim data overlaid with contours from the Achtermann & Lacy 
(1995) [Ne II] line map deconvolved with the maximum entropy method (their Figure 5). 



- 12 - 




